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Abstract 

Stationary perfect-fluid configurations of Einstein's theory of gravity are studied. It is 
assumed that the 4^-velocity of the fluid is parallel to the stationary Killing field, and also that 
the norm and the twist potential of the stationary Killing field are functionally independent. It 
has been pointed out earlier by one of us (I.R.) that for these perfect-fluid geometries some of 
the basic field equations are invariant under an SL{2JR) transformation. Here it is shown that 
this transformation can be applied to generate possibly new perfect-fluid solutions from existing 
known ones only for the case of barotropic equation of state p + = 0. In order to study the 
effect of this transformation, its application to known perfect-fluid solutions is presented. In 
this way, different previously known solutions could be written in a singe compact form. A new 
derivation of all Petrov type D stationary axisymmetric rigidly rotating perfect- fluid solutions 
with an equation of state p + 3p ~ constant is given in an appendix. 

PACS number: 04.20.Jb, 04.40.+C 
1. Introduction 

In spite of great efforts, no stationary axisymmetric perfect-fluid solution of Einstein's 
equations has been found which would be appropriate for describing, in the framework of 
general relativity, rapidly-rotating compact massive stars. In fact, only a small number of 
stationary axisymmetric perfect-fluid solutions without higher symmetries are known [1-8] (see 
also references therein). On the other hand, owing to the nonlinear character of Einstein's 
equations, only the knowledge of as much solutions as possible could provide us with a deeper 
insight into the methods that might help in obtaining astrophysically relevant solutions of 
Einstein's equations. Therefore, every new solution generating method is welcome in general 
relativity. 
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Here we would like to present and apply a technique by which one can generate new 
solutions from previously known ones. This technique is, in fact, a generalization of a method 
given by Geroch [9] by which out of any source-free solution of Einstein's equations possessing 
a non-null Killing vector field one could get a new one-parameter family of vacuum solutions. It 
has been recently shown by one of us [10] , that an analogous transformation can be successfully 
used for spacetimes possessing a non-null Killing field and certain kinds of matter fields. In 
particular, some of the basic field equations for perfect-fiuid matter sources having 4-velocity 
parallel to a timelike Killing field possess exactly the same type of invariance as the vacuum field 
equations [10,11] ^. A perfect fluid with 4-velocity parallel to a timelike Killing fleld is "rigid", 
i.e., it is expansion- and shear-free. Consequently, such a fluid seems to be far from being 
appropriate for astrophysical applications. However, as it was shown by Geroch and Lindblom 
[15], perfect fluids of this kind represent equilibrium configurations of dissipative relativistic 
fluids. Thereby, the study of these models is of obvious astrophysical interest. Note that this 
"rigidity" assumption introduces considerable simplifications of the basic field equations. 

In this paper we consider stationary perfect-fluid conflgurations with 4-velocity parallel to 
the stationary Killing fleld. The norm and the twist potential of the Killing field are assumed 
to be functionally independent. In section 2 we recall some of the basic notions and results of 
the projection (or 3-dimensional) formalism of general relativity developed for spacetimes with 
a non-null Killing field. Then, we present the field equations in geometrically preferred local 
coordinates. In section 3 a transformation by which one can generate new solutions from known 
ones will be recalled. In particular, we determine the conditions under which this transformation 
can be applied to the selected perfect-fluid spacetimes. It turns out that the original perfect 
fluid has to possess the barotropic equation of state p + 3p = which is invariant under the 
action of the transformation. Finally, in section 4, wc apply the developed technique to all 
existing suitable stationary perfect-fluid solutions. Accordingly, in section 4, first wc recall 
the most general family of known stationary axisymmctric perfect-fluid solutions for which the 
equation of state is p + 3P = 0, the flow is rigid in the above sense and the twist potential and 
the norm of the Killing field arc functionally independent. These solutions belong to the large 
family consisting of all Petrov type D stationary axisymmetric rigidly-rotating perfect-fluid 
metrics with the equation of state p + 3P = constant, obtained by Senovilla [Ij. Note that this 

^ The invariance properties of the basic field equations for electrically charged rigid perfect fluids were studied earlier by 
Kramer ct ol [12, 13]. Note, however, that their considerations were restricted to the StdtZC case exclusively. An application of 
their approach was the derivation of a regular static charged interior Reissner-Nordstrom solution from the interior Schwarzschild 
metric [14]. As far as we know, the results used in the present paper to analyze the invariance properties of the field equations 
are new and they were published first in Refs. [10,11]. 
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family contains as a subfamily each previously known stationary axisymmetric perfect-fluid 
geometry with the equation of state p + 3p = constant. Several other solutions were found 
later by Mars and Senovilla which are also of Petrov type D stationary perfect-fluid solutions 
- one of them possesses two commuting timelike Killing fields while the other is a stationary 
axisymmetric solution - and the equation of state is p -|- 3p = for both of them. In appendix 
1 these solutions are shown to be subclasses of the ones given by Senovilla. Finally, a new 
form of all Petrov type D stationary axisymmetric rigidly rotating metrics with an equation of 
state p + 3p = constant is presented in appendix 2. The solutions yielded by the application 
of the generalized Geroch transformation to the original Senovilla metrics can be related easily 
to these general but relatively compact and simple metrics. Moreover, this family of solutions 
contains on equal footing both of the original classes of the Senovilla solutions. 

2. Basic notions and the field equations 

Wc consider perfect-fluid spacctimcs, (M.gab), admitting a timelike Killing field, . The 
mass density, the pressure and the 4- velocity of the fiuid are denoted by p, p and u°- , respectively. 
We assume, furthermore, that the 4- velocity is parallel to the timelike Killing field which 
implies 

= {-v)-^i''. (2.1) 

Then, the flow is rigid in the sense that it is expansion- and shear-free. As mentioned above, 
this assumption on rigidity is compatible with the general strategy that one would like to find 
a faithful description of equilibrium states of dissipative relativistic fiuids. 

It is known that the entire geometrical content of a spacetime possessing a timelike Killing 
field can be represented on the 3-dimensional space of Killing orbits by a triple [v, iVa, hab) 
where v is the norm of the Killing field, 

V = rCa, (2.2) 

Wa is the twist of the Killing field, 

UJa = Cabcd^'^V''^'^, (2.3) 

(here Cabcd is the 4- volume element), while hab is a Riemannian metric defined as 

hab = -VQab + CaCfe- (2-4) 

For the selected stationary perfect-fluid configurations Ua can be expressed (at least locally) 
as the gradient of a function, u, called the twist potential 

Wa = DaW. (2.5) 
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Throughout this paper we shall assume that the above defined functions v and cu are 
functionally independent, i.e., 

D[„u;D^v ^ 0. (2.6) 

(The complementary case, when v and uj are functionally dependent, will be considered in a 
subsequent paper.) Whenever eq. (2.6) is satisfied, the basic field equations can be simplified 
by the introduction of geometrically preferred local coordinates, {x^,x'^,x^) [10,11]. Then we 
have that the functions v and uj and also - by virtue of the field equations and the equation 
of state which is of the form p — p[p) for the considered rigid perfect fluids the functions 
p and p depend only on the coordinates and x'^. In these local coordinates the set of 
equations describing a stationary perfect fluid with 4-velocity parallel to the stationary Killing 
field reduces to [10,11] 

2{Rab - ^^tiAB) = v~'^{dAvdBV + dAUJdBti>), (2.7) 

Rsp = Wnv-^phsp, (2.8) 

dAP+l{p + p)—=0, (2.9) 
2 V 

where Rat is the Ricci tensor associated with the three dimensional metric hab and the capital 
Latin indices take the values 1, 2 while the Greek ones take 1, 2, 3. As it follows from the results 
of Refs. [10,11], for the stationary perfect-fluid geometries under consideration it is sufiicient to 
solve the above (reduced) set of fleld equations, because the functional independence of v and 
CO ensures that any solution of this set satisfies the complete system of Einstein's equations. 

3. The transformation 

In this section first the symmetry properties of eq. (2.7) will be studied. In particular, it 
will be recalled that a one-parameter family of new solutions of eq. (2.7) can be associated with 
any given solution of this equation. Then, the restrictions on the applicability of the related 
transformation - imposed by the remaining fleld equations, (2.8) and (2.9) - will be determined. 

To start off, note that, just like for the vacuum case, the left hand side of eq. (2.7) depends 
merely on the 3-dimensional metric, hat, while the right-hand side is given in terms of the 
functions v and u). By exactly the same argument as used for the vacuum case in Ref. [9], 
it can be shown that by starting with a particular solution, (uo,wo), of eq. (2.7) associated 
with a fixed 3-dimensional metric, hab, (and thereby, a flxed set of functions Rab — ^^^ab) 
one can generate a one-parameter family of solutions, {vT,u)r), of eq. (2.7). More precisely, let 
(ucwo) be a particular solution; then, the full set of solutions of eq. (2.7) associated with a 



fixed S-metric, hat, can be given as 

Vo 



(cos T — u>o sin r)2 + Vq sin^ r ' 



(3.1) 



{uo cos T + sin r) (— wo sin r + cos r) — Vq sin t cos r 

'^■^ = ? '■ — 2~^"2 • (3.2) 

(cos r — Wo sin r)"" + sin r 

These formulas, with r as the only independent parameter, can be derived from the general 
form of an 5^(2,11) transformation by factoring out with respect to the trivial two-parameter 
gauge transformations [9] . Although this transformation may yield gauge inequivalent geome- 
tries, which can be shown to form a circle, its repeated applications do not generate new 
solutions since the transformation then merely yields a rotation of this circle ^. 

The most significant diflference between the set of vacuum field equations and the field 
equations for the case under consideration is that for the vacuum problem the appropriate 
form of eq. (2.7) is the only equation to be solved, while for the perfect-fluid case the basic 
field variables, beside eq. (2.7), have to satisfy eqs. (2.8) and (2.9), as well. These equations, 
however, impose non-trivial restrictions on the applicability of the transformation. 

To determine these restrictions, let us consider first eq. (2.8). It can be easily seen that if 
one fixes the 3-geometry, hab, the quantity pv~^ has to be left intact by the above transforma- 
tion. While V transforms according to eq. (3.1), p has to transform so that the variation of p 
compensates the variation of v. This implies that 

^ = ^, (3.3) 

Vr Vo 

has to hold, where po and vq denote the pressure and norm of the Killing field for the original 
perfect-fiuid configuration, while Pr and Vr are the corresponding functions for the transformed 
geometries. This equation represents one of the subsidiary conditions to be satisfied whenever 
one applies the above transformation. 

A further restriction is risen by the Euler-Lagrange equation (2.9). In fact, eqs. (3.1) 
and (3.3) can be used to determine the transformed mass density, Pr, and, thereby, the trans- 
formation law for the equation of state of the fluid. To see this, take the derivative of eq. 
(3.3) 

Vr [oaPt-Pt )='*^o [OaPo-Po )• (3.4) 

Vr Vo 

Then, using the Euler-Lagrange equation (2.9), one gets from (3.4) 

{pr + 3pr ) = (po + 3po) • (3.5) 
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The properties of the transformation are analyzed in detail in Ref. [9] for the vacuum case. Exactly the same analysis 
applies for the present case. 



Finally, substituting the right hand side of (3.1) for Vj- into (3.5) one obtains, 

(po + 3po)9awo = 

= {p-r + 3Pt-)| [(cost — Wo sin r)^ — sinr^i^oj^AVo + 2sinrwo(cosr — wq sinr)9^a;o|. (3.6) 

Since vq and wq are supposed to be functionally independent, the coefBcients of OaVo and 
Oai^o in eq. (3.6) must vanish separately. This, however, implies that both of the quantities 
Po + 3po and Pr + 3pr must vanish identically. Correspondingly, the above transformation can 
be applied only to those perfect-fluid configurations that have the barotropic equation of state 
p + 3p = 0. Furthermore, each of the new solutions yielded by the transformation has to possess 
this equation of state. 

It is important to emphasize that there is no further restriction on the applicability of 
the above transformation (see eqs. (3.1) and (3.2)), risen by eqs. (2.8) and (2.9). In fact, for 
these perfect-fluid geometries with equation of state p + 3p = the Euler-Lagrange equation 
can be integrated. It is easy to check that the general solution of eq. (2.9) for this equation of 
state is pv~^ = constant which, by virtue of eqs. (2.8) and (3.3), implies (in accordance with 
the fact that hab is kept fixed) that PrV^^ = PqVq^ = constant. Hence, for these perfect-fluid 
configurations the functional forms of Pr and Pr will not be explicitly presented. They can be 
obtained by simply multiplying Vr by constant factors. 

4. Getting solutions from existing ones 

In this section the above transformation is applied to known stationary perfect-fluid solu- 
tions of Einstein's equations with equation of state p + 3p = 0. The main steps of the procedure 
of getting new solutions are the following: Start with a stationary perfect-fluid solution with 
equation of state po + 3po = 0. Determine, first, the functions Vq and ojq as well as the fixed 
3-metric, hab, which are the input data for our procedure. Then, determine the functional form 
of Vr and Ur using eqs. (3.1) and (3.2). Finally, the 4-dimensional line element can be given, 
in virtue of eq. (16.22) of Ref. [13], by the integration of the relevant form of eq. (16.23) of 
Ref. [13]. 

The stationary perfect-fluid solutions which can be used as input for the above procedure 
are, in fact, rare. There is a large class containing Petrov type D stationary axisymmetric 
perfect-fluid solutions with the equation of state p + 3p = constant, given by Scnovilla [1], 
which includes as a special subcase the Wahlquist solution [13] and certain families of solutions 
given earlier by Kramer [5,6] and by Mars and Senovilla [3] 

It was not mentioned originally by the authors in Ref. [3] that two of the solutions given there represent, actually, a 
subclass of the previously published Senovilla metrics. 
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The general form of all Petrov type D stationary axisymmetric rigidly rotating perfect-fluid 
metrics with an equation of state p + 3p = constant can be written [1] as 

ds^ = vo{dt + Aod(pf + {V- W)[G-^dx^ + H'^dy^ + c^GH{G - Hy^dcp"^] (4.1) 



where the norm of the stationary Killing field, (^)", is 

H-G 
V -W 



= (4.2) 



furthermore, 

_ c{HV-GW) 

- H-G ' ^^-^^ 

and c is an arbitrary constant. Here the functions V = V{x), W = W{y), G = G{x) and 
H — H{y) satisfy eqs. (16) and (17) of Rcf. [1] (sec also eqs. (A2.1) and (A2.2) of the present 
paper). The Senovilla metrics are divided into two classes. The first one, class I with ^ 7^ 0, 
coincides for a special setting of the parameters [1] with the Wahlquist solution [13], while in 
class II with ^ = 0, one can also find metrics given previously by Kramer [5,6] and by Mars 
and Senovilla [3]. For the Senovilla metrics the functions V, W, G and H satisfy eqs. (18), 
(19a) and (19b) of Ref. [1] for class I and eqs. (18), (20a) and (20b) of Ref. [1] for class II. 
The general solutions of these equations for the case of p + 'ip ~ (which corresponds to the 
vanishing of the parameter a in these equations) can be given as follows: 
For class I: 

V{x) =mo + M sin(46a; + xq) (4.4) 
W{y) =mo + M cosh(46t/ + yo) (4.5) 
G{x) = so + 5sin(46x + xo) + Ni cos(46a; + xo) (4.6) 

H{y) = so + N2 sinh(46y + yo) + S cosh(46y + yo) (4.7) 

with M = ^Vm"^ - 'ib'^n, mo = .so = S = ^2^-^^^, where m, n, s and h are the 
constants used in the field equations of Ref. [1] and, finally, xq, j/q, A''i and N2 are constants of 
integration 



^ Note that the functional form of W( yj given by eq. (4.5) differs from the original expression of W(y) in Ref. [1]. After 
checking the functional form of W(yJ given by eq. (21) of Ref. [1] we have found that it does not satisfy eq. (18) of Ref. 
[1] unless the second term of the right-hand side of eq. (21) of Ref. [1] is divided by the constant Ci used there. Since the 
function y) depends on the explicit form of W(y) (see eq. (19b) of Ref. [1]), to have the correct solution given by eq. (4.7) 
for class I, one has to alter the function H^y) correspondingly. Furthermore, in the last term of the expression for G(x) for 
class I in Ref. [1] instead of the factor {Ahx -\- C4) there should be (-46a; + d) in order to have the correct solution of 
the basic field equations. 



For class II: 



m 



V{x) = = constant 



W{y) = L^e-^'>y + 



m 
262 



G{x) = — Z/3 cos(46a;) — Li sin(46a;) 



H{y) 



1662 



(4.8) 

(4.9) 

(4.10) 

(4.11) 



where Li,L2,L3 and L4 arc arbitrary constants. 

Using eqs. (4.2)-(4.7) the norm of the Kilhng field, vq, and the twist potential, ojq, for class 
I can be shown to be 



S + 



Ni cos(46x + xq) - N2 sinh(46y + j/o) 



sin(46x + Xq) — cosh(46t/ + t/o) 

1 N2 cos(46a; + xq) + Ni sinh(46?/ + yo) 
M sin(46x + a;o) — cosh(46y + yo) 

The same quantities for class II are 

vo = -^(^3 cos(46x) + J:4sin(46a;))e''*'^ - ^, 
Li Li 

Wo = ^e'^^y (-L3 sin(46a;) + L4 cos(46a;)) . 



(4.12) 
(4.13) 

(4.14) 
(4.15) 



Here the twist potentials have been determined from the relevant form of eq. (16.23) of Ref. 
[13]. 

The families of perfect-fluid configurations obtained from the Senovilla-metrics are given 

as 



dSr^ = Vridt + Ard(j)f 



1 



(4.16) 



G H 

where G and H are the original fimctions given by eqs. (4.6)-(4.7) and eqs. (4.10)-(4.11) for 
class I and class II, respectively. Moreover, the actual forms of Vr and lOt can be calculated by 
using eqs. (3.1) and (3.2) with vo and loo determined by eqs. (4.12) and (4.13) for class I and 
by eqs. (4.14) and (4.15) for class II. Finally, the function At = ^^(a;, y) is the solution of the 
system 

(4.17) 



OX oy 



dy ^ dx ^ 

Using expression (3.1) for Vr, it turns out that the general solution of eqs. (4.17) and (4.18) 
can be written as 



Ar = 



HWr - GVr 

H-G 



(4.19) 



where Vj- and Wj- are the transformed versions of the functions V{x) and W{y). It is also true, 
that Vr can be expressed as Vr = {H — G)/{Vr — Wr). The functions Vr and Wr are 



Vr{x) = /(r) + Cir sm{4bx + xq) + C2t cos{4bx + xq) (4.20) 

Wr{y) - fir) + sinh(46y + yo) + C4. cosHAby + yo) (4.21) 

where /(r) is an arbitrary function of t with the property /(O) = mo while the factors Ci,-, 
C2t, Cst and C4T are given by the following expressions. 
For class I: 

Ci^ = i- [sin2 t{S^ - Nf - iV|) + M2 cos2 r] , (4.22) 



2 sin. T 

C2r = — ^ [sin tNiS - cos tTVzM] , (4.23) 
2 sin T 

Cs^ = — ^ [sin TN2S + cos rA^iM] , (4.24) 
C4r = ^ [sin2 t{S^ + iVf + iV|) + cos" r] . (4.25) 



For class II: 



2 sin T 

Cir = (i2-^4SinT + iiL3COST), (4-26) 

Li 
2 sin 7" 

C2t = r {L2L3SmT - L1L4COST), (4.27) 

C3. = ^^^^i^sinV-Li, (4.28) 

C4. = -^^^ii^sinV + Li. (4.29) 

It is straightforward to check that for both classes I and II the identity = Ci^+C2r+C^r 
holds. Using this relationship, it can be seen that the transformed solutions Vr, Wr, G and 
H also satisfy the original field equations (16) and (17) of Rcf. [1] ^. Therefore, we can state 
that the transformed metrics represent Pctrov type D stationary axisymmctric rigidly-rotating 
perfect-fluid configurations, i.e. the solutions obtained by the application of the transformation 
arc, actually, not new. In particular, it can be shown that the metrics obtained from class I 
include (for special settings of the parameters) both of the original classes (^ 7^ and ^ = 0) 
of the Scnovilla metrics, while the metrics obtained from class II comprise a subclass of this 
family. Note that although the transformation docs not yield new solutions from the Scnovilla 
metrics, it is not a gauge transformation. This follows from the fact that it maps both class I 
and class II onto a strictly larger subclass of the entire family of solutions. 



See also the general form of the functions V, W,G and H given in appendix 2. 
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5. Conclusions 



Stationary perfect-fluid configurations with 4- velocity parallel to the stationary Killing field 
were considered for which the norm and the twist potential of the Killing field arc functionally 
independent. It was shown that by the use of 'effective' SL{2,1R) transformations new perfect- 
fluid solutions can be obtained from known ones whenever the equation of state is p + 3p = 0. 
By applying the relevant method to known geometries with this equation of state we obtained 
stationary perfect-fluid solutions of Einstein's equations which contain in a compact form all 
the configurations we used as input for our procedure. Although these solutions are not gauge 
related to the ones we started with, there are no new perfect-fluid configurations among them, 
since all special cases were given previously. 

Note, finally, that the perfect-fluid solutions to which our procedure has been applied here 
possess at least two commuting Killing flelds. However, we would like to emphasize that our 
method postulates merely the existence of a single timelike Killing field. 

Appendix 1 

Here it is shown that the two families of metrics given by formulas (27) and (42) of Ref. 
[3] are, in fact, subclasses of class II of Senovilla metrics. 

The case of a Mars-Senovilla metric with two timelike Killing fields 

A stationary rigid perfect-fluid solution with equation of state p + 3p = was found by 
Mars and Senovilla [3] . This spacetime is also of Petrov type D and its line element is 

^ 2 ,.2 . jrr. a .2 a'^y'^dx^ ^ /^i in 

ds = ^r^dt -[aydT+—dt) +— — + dy , (Al.l 

a^y^ a/3' ay ' - ^^{x^ - 2f3^) 

where a > 0, (3 > 0, 5 and a are arbitrary constants. This solution possesses two commuting 
Killing vector fields, (^)" and (^)°, associated with the coordinates t and T, respectively. It 
is an interesting feature of this solution that both of the commuting Killing fields are timelike. 
The 4- velocity of the fluid is parallel to the Killing field (^)"- 

To show that this metric is isometric to class II of Senovilla metrics change the variables 

as 

lex / — / — 
dt ^ J^dt , dx^ VGdx , dy VHdy , d(j) dT {Al.2) 

and set the functions in the Senovilla metrics to 



G{x) = Q^{x), V{x) = 0, H{y) = 0^ = constant, W{y) = -e'^"^, (A1.3) 

10 



where a and /3 are constants of the Mars-Senovilla metric. Substituting the above expression 
for G{x) in eq. (20a) of Ref. [1] and changing the constants of Ref. [1] to 

we get the corresponding equation of Ref. [3] for Q{x) as given at the bottom of page 3063 in 
Ref. [3]. Prom this point, the same way as in Ref. [3], one obtains the metric in the form given 
by eq. (42) of Ref. [3]. 

Finally, note that the "axial" Killing vector of the Senovilla metrics, with the settings 
(A1.2)-(A1.4), becomes, actually, a timelike Killing vector. 

The case of a stationary axisymmetric Mars-Senovilla metric 

Another stationary axisymmetric Petrov type D perfect-fluid solution with equation of 
state p + 3p = was also found by Mars and Senovilla (see case I of Ref. [3]). Originally, 
this solution was thought of and described as a differentially rotating perfect fluid. However, 
the authors informed us, that, although the line element is correct (see eq. (27) of Ref. [3]), 
this stationary perfect-fluid solution is, in fact, rigidly rotating with 4-velocity parallel to the 
stationary Killing field. 

The line element of this spacetime is 

where (3, u, a and Lq are arbitrary constants and e is a sign [3]. Here is a function of the 
variable y satisfying 

H = ea^H, (Al.6) 

where the dot denotes derivative with respect to y. 

To show that this metric can be obtained from class 11 of Senovilla metrics, first, change 
the variables as 

dt , dx \/Gdx , dy ^/H^dy (^1-7) 
VP 

in the Senovilla metrics where Hg{y) is the function given by eq. (4.11). Furthermore, set 

Gix)=u;'L^{x), V{x) = 0, /f^(y) = |JM, W{y) = -H-^{y), {A1.8) 

where (3 and cu are constants of the Mars-Senovilla metric. Then, setting the constants of the 
Senovilla metrics as 

/^ = 4ea^ 6=^, c=-, {A1.9) 
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one can see that the functions G{x), V{x), Hg{y) and W{y) satisfy eqs. (18) and (20a-b) of 
Ref. [1] whenever the function L{x) satisfies eq. (25) of Ref. [3] and H{y) satisfies the equation 
given bellow eq. (27) of Ref. [3]. 

Appendix 2 

Here we are going to reflect briefly on a method yielding solutions of the field equations 
(16) and (17) of Ref. [1]. This way one gets the general form of all stationary axisymmetric 
Petrov type D rigidly rotating perfect-fluid metrics with an equation of state p + 3p= constant 
possessing a simpler functional form than the original solutions given in Ref. [1] . 

The metric is given by eq. (4.1) and the basic field equations for the functions V = V{x), 
W = W{y), G = G{x) and H = H{y) (see eqs. (16) and (17) of Ref. [1]) are 

Y>2 ^ ^ 2{V - W)[8b^iV -W) + V"], {A2.1) 

12a{V -W) = l%h\G -H) + - + G", {A2.2) 

where the prime denotes derivation with respect to the argument. The general solutions of 
these equations can be obtained as follows. 

Taking the derivative of eq. (A2.1) with respect to the variable x and using the fact that 
for a non-singular metric V — W we get 

{V" + 166V) = 0. (A2.3) 

The general solution of this equation is 

V{x) = Vo + Ci sm{Ahx) + C2 cos(46a;), (A2.4) 

where Vq, Ci and C2 are constants. 

By a similar procedure, calculating the derivative of eq. (A2.1) with respect to y, we get 

W'{W"' - mb'^W) = (A2.5) 

which holds either if 

W = (A2.6) 

or 

W" - Wb'^W' = 0. {A2.7) 
The general solution of eq. (A2.7) for the case when W 7^ is 



W{y) = Wo + C3 sinh(46y) + C4 cosh(46j/). 
12 



{A2.8) 



The substitution of expressions (A2.4) and (A2.8) into eq.(A2.1) yields the following re- 
strictions on the constants 

Vo = Wo =: Co and C| = Cl + C| + C|. (A2.9) 

Substituting W = and eq. (A2.4) into eq. (A2.1) we have 



W{y) = Co ± \Jcl + Cl = const . {A2.10) 

To get the solutions to eq. (A2.2), first note that by using 

V" + 166V = mb^Co {A2.11) 

eq. (A2.2) can be recast into the form 

12a{V -W) = m^{G - H){Co -W)- V'G' - W'H' + G"{V - W). {A2.12) 

Taking the derivative of eq. (A2.12) with respect to x, using (A2.11) and the fact that 
V — W ^0 we arrive at 

G'" + IQb^G' = 24aV'. (A2.13) 
This equation has the general solution 

G{x) =Go + Ki sm{Abx) + K2 coa{Abx) - , (A2.14) 

where Gq, Ki and K2 are constants of integration. 

In order to determine the function H{y) we take the derivative of eq. (A2.12) with respect 
to y and use 

w" - m'^w = wb'^Co. {A2.m) 

This yields 

24a{V - W)W' = W'{G" + Wb'^G + H" - 16b'^H). {A2.16) 
First, consider the case when W' ^ 0. Then from eq. (A2.16) it follows that 

H'" - m'^H' = -24aW' {A2.17) 

the general solution of which is 

H{y) = Ho + Ks sinh(46y) + cosh(46j/) - ^yW, (^2.18) 

where Hq, K3 and K4 are constants of integration. 
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By the substitution of expressions (A2.14) and (A2.18) for G{x) and H{y) into eq. (A2.15) 
one can see that the constants of integration have to satisfy the following conditions 



Ho = Go=:Ko and K4C4 = K^Ci + K2C2 + K3C3. (A2.19) 
Finally, whenever W = eq. (A2.16) is satisfied identically. Then the substitution of 



W = Co± ^/C'i + cl into eq. (A2.12) yields 

C K + C K 

H{y) =Ko± ^ = constant . {A2.20) 

Our solutions can be related to the original Senovilla ones by setting the constants m, n, 
s and h, introduced in the field equations (18)-(20) of Ref. [1], in terms of our constants of 
integration as follows: 

m = 2b^Co, s = h'^Ko - 2aCo, n = b^iC^ - - C|), 

3 (^2.21) 

h = a{C^ - Cl - Cl) - -b^CoKo - CiKi - C2K2). 

Finally, note that the class II Senovilla metrics can be recovered from our solutions by 
choosing the constants Ci = C2 = 0. 

Acknowledgments 

One of us (I.R.) wishes to thank J. M. M. Senovilla for useful discussions and also for 
drawing into our attention a number of solutions which could be used in our formalism. Thanks 
are also due to A. Sebestyen for helpful discussions. A great number of calculations has been 
carried out using the algebraic computer programme MAPLE V and also by making use of 
GRTensor II [18]. This research was supported in parts by the OTKA grants F14196 and 
T016246. 



14 



References 

J. M. M. Senovilla: Phys. Lett. A 123 number 5, 211-214 (1987) 

J.M. M. Senovilla: Class. Quant. Grav. 4, L115-L119 (1987) 

M. Mars and J. M. M. Senovilla: Class. Quantum Grav. 11 3049-3068 (1994) 

S.Bonanos and E. Kyriakopoulos: Class. Quantum Grav. 11 L23-L28 (1994) 

D.Kramer: Class. Quantum Grav. 1 L3-L7 (1984) 

D. Kramer: Class. Quantum Grav. 2 L135-L139 (1985) 

E. Kyriakopoulos: Gen. Rel. Grav. 20, 427-436, (1988) 

F. J. Chinea and L. M. Gonzales-Romero: Class. Quant. Grav. 1271-1301 (1992) 
R. Geroch: J. Math. Phys. 12, 918-924 (1971) 

I. Racz: ESI Preprint Series No. 168/94 (1994) 

I. Racz: in Inhomogeneous Cosmological Models, eds. A. Molina and J. M. M. Senovilla (World 
Scientific, Singapore, 1995) 

D. Kramer, G. Neugebauer and H. Stephani: Fortschr. Physik, 20 1-75 (1972) 
D. Kramer, H. Stephani, M. MacCallum and E. Herlt: Exact solutions of Einstein's field equa- 
tions (Cambridge; Cambridge University Press, 1980) 

D. Kramer, G. Neugebauer: Annalen Physik 27 129-135 (1971) 
R. Geroch and L. Lindblom: Ann. Phys. 207, 394-416 (1991) 
C. M. Cosgrove: J. Phys. A 11, 2389-2404 (1978) 

E. D. Fackerell and R. P Kerr: Gen. Rel. Grav. 23, 861-876 (1991) 

P. Musgrave, D. Pollney and K. Lake: GRTensor II, (Queens University, Kingston, Ontario, 
Canada 1994) 



15 



